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Abstract 

A new mathematical notation is proposed for the iteration of functions. 
It facilitates the application of the iteration of functions in mathematical 
and logical expressions, definitions of sets, and formulations of algorithms. 
Illustrations of the notation include definitions of constant points, periodic 
points, a filled-in Julia set, the Mandelbrot set, iterations of a logistic map, 
the double-approximating procedure for solving the Lorenz equations, a 
description of a financial time series, and reordering nonnegative integers 
useful for the investigation of the Collatz's (3x+l)/2 convergence problem. 
The terms iteral and iteral of function are suggested to name the new 
denomination. 



1 Introduction 

A mathematician or a natural science researcher routinely writing in equations 
f{x), nl, ^2, FJ, and J relies on clarity, expressiveness, and convenience of math- 
ematical notation filtered out by centuries of work, recognition, and acceptance. 
Creating a theory in many respects is developing a language eventually accom- 
panied by a system of signs and symbols. A historical excursus ^ pp. 441 - 443] 
and (3) shows that the road of many habitual mathematical symbols was bumpy 
with periods of oblivion, sometimes rediscovery, and final failure or success. A 
lack of a well understood notation creates a barrier for a wider application of 
a subject. In attempt to reuse the iteration of functions for the analysis of the 
financial high frequency time series the author did not find a notation suitably 
expressing this concept for writing formulas and equations. 

Latin iterum means again. It is the root for the word iteration. Iteration is 
a fundamental notion in computer science, where it denotes repetition of steps, 
the repetition of a sequence of instructions in a computer program or algorithm. 
Accordingly, all computer languages have the means to express iterations. For 
instance, in C++ [27, pp. 136 - 137] these are for, while, and do-while 
iteration statements. The Unified Modeling Language abstracts the iterative 
expression and places the iteration above concrete programming languages |191 
pp. 318 - 319]. 

Dynamical systems, theory of chaos, and fractals depend on iterations. Ben- 
oit Mandelbrot comments [161 p. 260, p. 268]: 

The motivation for iteration and many of its important properties come 
straight from physics. In iteration, physicists see a simplified view of the state 
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of a dynamical system and its evolution in time. . . . The theory of iteration of 
rational functions goes back to the mid-nineteenth century and perhaps even to 
Abel. 

Marek Kuczma et al. [11, p. 13] emphasize the relationship between itera- 
tions and functional equations: 

We feel that any attempt to divorce iteration from functional equation in- 
vestigations would be an extremely, indeed totally, fruitless task. 

This article proposes a new notation and term for iteration of functions and 
illustrates their usefulness. 



2 Existing Notations Implicitly Using Iterations 

The sign ^2 for summation was introduced by Leonard Euler in Institutiones 
calculi differentials, St. Petersburg, 1755, Cap. I, $ 26, p. 27, who was saying, 
summam indicabimus signo ^ [3J item 438, p. 61]. In its variations the notation 
assigns initial values to the involved quantities, defines the number of iterations, 
binds each iteration with an arithmetic addition of the current term to a previous 
accumulated result, and assigns a value to the operation. It is unambiguous. 
What is going on in the following left-hand sides is intuitively well understood 
and explained by the middle expressions 



n(n + 1) 



Al_ = l 1_ 1_ f_ 

^ i 2 ~ 1 + 2 2 + 3 2 + " ' ~ 6 ' 

i=i 

The expressions assume values, right-hand sides, and can become parts of other 
valid mathematical expressions. Presence of the counting index in the expression 
under the sum is not essential. Occasionally, the number of iterations is equal 
to zero. Then, a default, often zero, value is set for the sum. 

Another successful sign for n-factorial of natural numbers denoting "n(n — 
l)(n— 2)....3.2.1" was introduced by Christian Kramp of Strasburg in his Elemens 
d'arithmetique universelle, Colonge, 1808 [3, item 448, p. 72]. By convention 
0! = 1. The finite number of iterations is bound to the multiplication of factors 
incrementing by one and starting from one. The left-hand side below is more 
convenient than the middle or right-hand side 



31 = 1x2x3 = ^1x^1x^1. 



i=i 



It can shorten descriptions: 29! vs. 8841761993739701954543616000000. 

The factors in n\ are fixed for natural n. The signs for the product of n 
arbitrary factors were introduced by Thomas Jarrett, who extensively studied 
algebraic notations and published results in 1831 [3J item 448, pp. 69 - 71]. His 
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simplest proposal commented as "the product of n factors, of which the mth is 
a m " was P,"a m . The n was exactly above P. Today we write 



In 1812 Carl Friedrich Gauss introduced Y[ m the designation of continued in- 
finite products pO item 451, p. 78]. The symbols Hz were used by Gauss in 
Werke, Vol. Ill, Gottingen, 1866, p. 146 for expressing the gamma function 
item 650, p. 272]. Camill Jordan applies the notation l\ r (c 7Tpr - K n ) e for 
the multiplication of factors resulting from different values of r published in 
1870. George Chrystal represents the product (1 + iti)(l + 112).... (1 + u n ) by 
Y[ n 1 + u n printed in 1889 [21 item 451, p. 79]. The home page of the Clay 
Mathematics Institute displays a photocopy of the first page of the hand writ- 
ten manuscript of Bernhard Riemann's paper Ueber die Anzahl der Primzahlen 
unter einer gegebenen Grosse dated by 1859, November and first time published 
in Monatsberichte der Berliner Akademie in the same month [4|. The uppercase 
cursive Greek letters II and £ used without indexes build the starting Riemann's 
point, the identity proved by Euler and known as the Euler product formula 



Here p ranges over all primes 2, 3, 5, 7, ... and n ranges over all positive 
integers 1, 2, 3, 4, what is frequently symbolized by placing p and n under 
the product and sum signs. The right-hand side contains infinite number of 
terms and diverges for s = 1. Since the left-hand side involves unique primes 
without repetitions it would converge for any finite number of factors. This 
implies infinite number of primes. Jay Goldman comments [51 p. 36]: 

This may seem a rather difficult way of proving a theorem that Euclid had 
already proved by very elementary reasoning, but in the 19th century Dirichlet 
vastly extended the ideas used here and they form the basis for the most important 
applications of analysis in modern number theory. 

For us it means that Riemann similarly to Gauss applies the product sign for 
the infinite number of factors. In [2J, [3], and |33| the author did not find, when 
and by whom the symbol Y[ was used first time for the product of n factors. 
However, the popular mathematical handbook of the second half of the 20th 
century lists it among the first mathematical denominations in introductory 
sections [TO]. The product with zero iterations is assigned the default value one. 

A well known notation combining the limit with the sum of terms is J. 
Most frequently it implies infinite number of iterations. Jacob Bernoulli was 
the first, who used the word integral in print. Johann Bernoulli claimed for 
himself the introduction of the word. The principal symbol of the integral 
calculus has been suggested by Gottfried Leibniz as the long letter J reminding 
about calculus summatorius, the name favored by Leibniz over the name calculus 
integralis and the capital letter / favored by Johann [3, item 539, pp. 181 - 182]. 



n 
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The sign J is successfully applied to denote deterministic Riemann, Stieltjes, 
Lebesgue, Riemann-Stieltjes, Riemann-Lebesgue [10], and stochastic Ito, and 
Stratonovich integrals [T51 p. 48, p. 106]. Frequently, in abstract formulations 
under the integral sign we find useful for us combination of characters f{x). 

The first appearance of / for "function" together with parentheses is dated 
by 1734 and reported for Euler's Comment. Petropol. ad annos 1734-1735, 
Vol. VII, 1840, p. 186. This is cited by a secondary reference of Johannes 
Tropfke in [3J item 643, p. 268] as the Euler's words: "Si /(§ + c) denotet 
functionem quamcunque ipsius - + c". Let us list the following common and 
useful properties of the discussed notations 

• The symbols X^^ML / i m ply finite or infinite iterations. 

• The number of iterations and initial quantities are defined. 

• The f(x) under J is visible and can be substituted by a concrete function. 

• The denominations assume values and can be parts of other expressions. 

• The notations are unambiguous and can be mixed. 

What prevents applying the symbols based on iterations for defining the 
iteration of functions is binding ^2 ancl / to addition and Y[ an d n\ to multi- 
plication. Let us review existing notations for the iteration of functions. 

3 Existing Notations for Iterations of Functions 

Articles and monographs are written on functional equations, dynamical sys- 
tems, chaos, fractals, and programming languages relying on iterations. The 
operation of composition is the basis of the key definitions. According to |11) . 
this is the only inner operation which can be defined in the family of all self- 
mappings T(X) of a given set X. All also includes such an element / £ F{X), 
denoted idx, that f(x) =iVi£ J. Since the composition o is associative and 
/ o idx = idx / V/ £ ^(X), the system (F(X),o) forms a semigroup with 
identity idx, monoid. In general, it is noncommutative. With these properties 
the iterates are defined as the powers /" of an element / £ T{X), where n is a 
positive or zero integer 

neN ,/°-id x ,/" +1 = /o/". 

While two arbitrary elements of the family do not have to commute with respect 
to the composition, the iterates of a function / <E F(X) do 

n,meN ,r°/ m = / m °/ n = f n+m . 

Here we are switching from the terms element and self-mapping to function 
using Euler's notation f(x). The composition of two functions f{x) and g{x) is 
expressed as (/ o g) = f{g(xj). 
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The x is a predecessor of y under / and y is a successor of x under /, if 
there exists an n € No such that y — f n (x). The ordering between predecessor 
and successor is partial. It is reflexive, and transitive but not antisymmetric. 
For example, if x 7^ 1 and f(x) = -, then b = f(a),a = f(b), and a^b. The x 
and y are equivalent, if they have a common successor under /: f n (x) = f m (y). 
This is an equivalence relation, where reflexivity and symmetry are obvious 
and transitivity follows from f m+k {x) = f m+l (y) = f l+m (y) = f l+n {z). The 
equivalence classes under / are called orbits [28, pp. 13 - 18]. This relation 
is fundamental for the theory of iterations. Gyorgy Targonski [28} p. 14, p. 
283] describes that it was introduced in a short remark made by Kazimierz 
Kuratowski in the text of a paper written in 1924 by Ralph Tambe-Lyche. It 
was independently rediscovered 17 years later and printed as a 14 lines abstract 
in 1941 by Gordon Thomas Whyburn, who coined the terms orbit and orbit 
decompositions and began systematic study of the new equivalence classes [ 1 1 L 
p. 14]. Possibility to select x, keep it intact, and change the power n e N 
sequentially creates the set {f n (x)} called the iteration sequence. Targonski 
|28j suggests to name it splinter of x borrowing the term from the theory of 
recursive functions, where it was introduced by Joseph Ullian [31J. 

The book [28] is considered in |11[ p. 14] as the first monograph completely 
devoted to the iteration theory. Systematically applying f n (x) and o it success- 
fully develops the theory proving notation's sufficiency. In practical applications 
the set X usually has an algebraic structure, where frequently it is needed to 
distinguish the iterates and powers with respect to multiplication. Anticipat- 
ing such a complication the authors of [TT| propose f 2 (x) for the composition 
f[f{x)] and {f(x)] 2 for the power f(x)f(x). Something is needed. Otherwise, 
what is the value of cos 2 x + sin 2 a; in a mixed environment? It is unity for 
powers and something else, if at least one term assumes composition. Following 
the proposal the mixed expression cos 2 x + [sin a;] 2 would be confusing because 
the parentheses are general grouping symbols. What have others applied in the 
historical order of publications? 

Olexandr (Alexander) Sharkovsky deriving properties of a continuous map- 
ping T of the line into itself denotes iterates T k a [22] • He switches to lower and 
upper indexes for an infinite sequence of iterated values of x, {T^x}°^ [23j . 

Tien-Yien Li and James Yorke |13] write x n +\ = F(x n ) and expand it for 
the logistic equation x n+ i = rx n [l — x n /K}. One recursive step is represented 
well. However, to use it in expressions, where initial values are assumed and the 
number of iterations is arbitrary is problematic. 

Mitchell Feigenbaum contributing to nonlinear dynamics applies the same 
x n+ i = f(x n ) for recursion equations [7J. The function f(x) can be easily 
replaced by a concrete one. However, substituting either the left-hand or right- 
hand side into other expressions either can be interpreted as vector coordinates 
or loses the number of iterations, when it is greater than one. 

Robert Devaney presents several programs for visualization of notions of 
chaos, fractals, and dynamical systems. Computer science demands iterations. 
Devaney exploits two step definitions like: for x > 1 and T(x) — x 2 ,T n (x) — > 
00 as n -> 00 I p. 12]. For instance, T 3 (1.2) gets the value 4.29981696. 



It is not clear how to replace T by x 2 . A reader can either be puzzled by 
(a; 2 ) 3 (1.2) or give the wrong answer 2.985984 for (1.2 2 ) 3 . Doing the composition 
in two steps is similar to writing P{x) = a + bx + cx 2 and J P[x)dx instead of 
J (a + bx + cx 2 )dx. 

Donald Knuth working with iteration of series favors x n 4— f(x n —i) for an 
iterative process and x n — /(/(.. ./(xo)...)) for the n-fold composition of a given 
function / [9, p. 530]. He chooses square parentheses f^(x) = (x)), 
reducing ambiguity with powers of multiplication. Sequential values Xji can 
be misinterpreted as vector coordinates. Replacing f^{x) by an actual func- 
tion can be misused. For the polynomial from the previous paragraph we get 
pN(x) = (a + bx + cx 2 )^ what is not intuitive. 

Mandelbrot considers the map z — > Az(l — z) and iterated maps z n = 
f(f(....f(zo))) [16, p. 38]. The last right-hand side is good for illustrations 
of iterates and not suitable for expressions. It is not concrete for referencing the 
number of iterations. Using the concrete left-hand side z n would be ambiguous 
in cases involving vector coordinates. 

We see that there is no common notation for the iteration of functions. 
Existing ones address theoretical needs and become ambiguous or less intuitive 
in mixed application environments. 



4 Proposal of Iteral 

When a mathematical notion is getting fundamental and other disciplines are 
demanding it, it may be a time for a common notation. The iteration of func- 
tions deserves it. The following are the simplest proposals 

CJVMRDT ITERATIONS PERFORMED f ( \ 
& I MDU% im VALUE J \ x l or 

CVA/mnT ITERATIONS PERFORMED f /\ 

a i ma(jij x=imTIAL value J \ x ) 

Here, SYMBOL should correspond to the meaning of iteration, INITIAL 
VALUE is an element of the domain set of the function f(x), ITERATIONS 
PERFORMED is the number of iterations completed so far and started from 0. 

If ITERATIONS PERFORMED = 0, then the value of the expression = 
INITIAL VALUE. A value of the function f(x) after a previous iteration, or 
INITIAL VALUE after the beginning, is a new argument of the function for the 
next iteration. Iterations may continue up to ITERATIONS PERFORMED 
only if intermediate values are members of the domain set of the function f(x). 
The second version explicitly indicates the variable x, where and expression for 
f{x) makes it unclear. This notation combines the following useful properties: 

• The function is visible and can be replaced with a concrete one. 

• The initial value is visible and can be replaced with a concrete expression. 

• The number of iterations is unambiguously specified. 

• The expression has a value and can be a part of other expressions. 
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• Combination of structural parts of ^ and / f{x)dx is intuitive. 

Greek letters are frequent candidates for symbols. To be unambiguous, they 
should differ from English letters. By this reason, Greek upper case iota is not 
suitable. Invention of a new symbol like J is challenging. The last reminds 
about summatorius and our should remind about iterum and iteration. Latin I 
has distinguishing dots but it is not as exact as I. Let us create a mirror image of 
N — > IL This makes it the first Cyrillic letter in Russian Hmepavflifi, translation 
of English Iteration, and a reminder. Thus, 

Wj{x) or VT x=v f{x) 

are the new notations. The author suggests to name this operation the iteral of 
function or simply iteral. This correlates with the formation rules of the words 
integral, differential, factorial, and fractal. Due to the Latin root, iteral can be 
naturally consumed by other languages. 

5 Applications of Iterals 

Below (. . .) and {. . .} denote a sequence and a set. Z, N, R, C are sets of integers, 
positive integers, real, and complex numbers. No = {0, N}. (No) = (0, 1,2, ...) 
is the infinite sequence, (n € N ) = (0, 1,2, ...,n) is a finite sequence, chain, of 
nonnegative integers. The set builder notation {x : P(x)} means the set of all x 
for which P(x) is true. The universal quantification of proposition V "for all", 
the existential quantification of proposition 3 "there exists", e "member of", A 
"and", V "or" are standard. 

Iterals have values. To get flexibility: H"x 2 = 2, H^x 2 = 4, H 2 x 2 = 16, 
Hfx 2 -> oo, YLfx 2 = 1, Hg° 5 x 2 -> 0. For x = f ,cos 2 x+ H 2 sin x = \ +sin^. 
Here are more iterals 

Tin jjn TT°° 1 Tin 1 9 

Yi Q (x + 1) = n, VV^ax = a n , Yi 1 = lim YY 1 = = 

° V 1 X ~ X 1 X+l n^oo 1 x+l 1 + V5 

The last is reciprocal of the golden ratio. Indeed, V(n £ N ) the iterals form 
the chain ("p^ifjfifjYjjj • ■ •)> where each member is the ratio of a previous 
Fibonacci's number to the current one. Using Targonski's and Ullian's term 
proposal, the splinter of 1 for the function is the set {H™^-} (n e No). The 
splinter is an infinite sequence, acyclic splinter, if Vm,n G N m ^ n H™/(x) 7^ 
Wff(x). The splinter is cyclic, if 3m, n e N to < n H" l /( x ) = H".f(x). 

Iterals form algebraic expressions and functions. Since the initial 
value of the iteral substitutes the argument of the function under the iteral 
sign the notation creates algebraic expressions or functions, if the initial value 
contains parameters and/or the function argument. The composition of func- 
tions must be valid. Iterals can nest. Here are a few examples used later for 
reordering N . 

Vn G No W v p = v = id*, H" p p = 2p, Yl n 2p+l p = 2p + 1, 
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Vfc G No H 2p+1 (2p + 1) = 2 k+1 p + 2 k+1 - 1, 

The last is proved by induction. For k — it holds: H^ p+1 (2p + 1) = 2p + 1 = 
2 0+1 p + 2 0+1 - 1. If it holds for k, then for k + 1 the left-hand side is equal to 

M*£ 1 (2p+ 1) = M2 fc+lp+2fc+1 _ 1 (2p+ 1) = 2 fe+2 p + 2 k + 2 - 1, 

and the right hand side is also equal to 2C £ + 1 )+ 1 p+2( fe + 1 )+ 1 -l = 2 k+2 p+2 k+2 -l. 

The proved iterals describe unfixed subsequences with the names OkO defined 
in the section on subdivision of No. The index k is the number of prefixing extra 
characters O, meaning odd, in the name of odd numbers O. For k = the name 
is just O, and the subsequence, given by the formula 2p + 1, represents all 
positive odd numbers with positions p = 0, 1, ... within this subsequence. The 
following will describe the fixed subsequences EOkO of the odd numbers 

Vi 1 2p y[ k 2p+1 (2p+l) = 2 k + 2 p + 2 k + 1 -l, k = 0,l, ... . 

The -E associates with even. The fixed subsequence of odd numbers with the 
shortest name EO for k = is defined by Ap + 1. Geometrically, the lines 
2 fc+2 p + 2 fc+1 — 1 of different orders /c do not intersect for p > 0. 

In general, nesting or composing iterals for the same function with common 
initial values is not the same as increasing the number of iterations 

M 2p+1 H2 P+I 2p - H 2p+1 (4p + 2) = 8p + 6, 

H 2p+1 2p= Hi p+2 2p = 8p + 4. 

This is because going from iteration to iteration the iteral's function remains 
intact and the initial value changes 

K +1 m - W\\ f , /;■<!. 

While nested iterals keep the common initial value intact and likely change the 
function on each composition step. 

We need more examples for the section on subdivision of N 

Vm e N H™ 2p = 2 m+1 p. 

For to = it holds by the definition of iteral H^ p 2p = 2p and 2 0+1 p = 2p. The 
induction step from m — >• m+ 1 gives the same result for the left- and right-hand 
sides: H™ +1 2p = ¥Ll m+lp 2p = 2 m+2 p and 2 ( - m+1 '>+ 1 p = 2 m+2 p. 

Vm, Z G No (tI 2p+1 • • • M 2p+1 ^ W 2p 2p - 2 m+l+1 p + 2 m+1 {2 1 - 1), 
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The index I is the number of composing iterals with the common initial value 
2p + 1 . Due to the previous proof the left-hand side is equal to 

^-* A 2p+l 1A -2p+l J L P- 

For I = the equality holds. In this case the iterals in the left-hand side vanish 
leaving 2 m+1 p = 2 m+0+1 p + 2 m+1 (2° - 1) = 2 m+1 p + = 2 m+1 p. The equality 
holds for / = 1 

Vll p+1 2 m+1 P = 2 rn+2 p + 2 rn+1 and 2 m+1+1 p + 2 m+1 (2 1 - 1) = 2 m+2 p + 2 m+1 
On the induction step I I +1 the left-hand side gives 

(Vl\ p+1 ■ ■ ■ Hl p+1 ) 2"^p = y[l p+1 (2^p + 2™+ 1 (2< 1)) = 
V / i+i 

= 2 m+l+2 p + 2 m+l+1 + 2 m+1 {2 1 - 1) = 2 m+l+2 p + 2 m+1 (2 l + 2 l - 1) = 

= 2 m+l+2 p + 2 m+1 (2 l+1 - 1). 

The same result is for the right-hand side 2 m +('+ 1 )+ 1 p + 2" l + 1 (2('+ 1 ' - 1) = 
2 m+i+2 p + 2 m+i( 2 ;+i _ ^_ Then; the foUowing result is simple 

Vm,l e N Wl p (W\ p+1 ■ ■ ■ H£2P = 2™+ l + 2 p + 2^(2' 1), 

For I > 1 the last iterals correspond to fixed subsequences of even numbers 
defined later and named EOiE m E. The shortest name EOE for I = 1 and 
m = corresponds to the sequence 8p + 2. 

Similar to tables of sums, products, and integrals it would be useful to have 
tables of iterals and add iterals into symbolic processing computer systems. 

Iterals work in definitions and statements. Fixed point v: f(v) = v. 
Using iterals: ¥l\f{x) = v. The original is shorter. However, the following 
property of a fixed point is expressed unambiguously with iterals: Vn G N 
H™/0) = v. Examples: Vn e N Vl^x 2 = 1 and H"x 2 = 0. 

Periodic point v of the order m: 3k, m 6 N Vfc < m H" l /(x) = v A H^f(x) ^ 
v. Example: Vw € R A v ^ H^i = v A 7^ u. 

j4 /i//ed in JitZza set. Biographical facts of Pierre Fatou and Gaston Julia 
and historical comments on iterations of rational functions, particularly, the 
complex function f(z) = z 2 + c, can be found in |16l pp. 268 - 275]. The iteral 
definition of the set is: z € C, J c = {v G C : \H^(z 2 + c)| < 00}. 

The Mandelbrot set. z € C, M = {c € C : IH^z 2 + c)| < 00}. The 
iteral definitions of the last two sets give straight forward directions for writing 
computer programs depicting the points of the complex plane. CH — h provides 
for the class complex making arithmetic and functional programming involving 
the complex numbers natural |27l pp. 267 - 274]. 

Iterations of logistic map. Two Belgian mathematicians Lamberte Adoplhe 
Jacques Quetelet and Pierre Frangois Verhulst contributed to modeling of the 
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growth of a biological population using an ordinary differential equation with 
separable variables of the population size and time. A solution of it is an S- 
shaped function. Ka-Kit Tung [29, p. 100] writes: 

Verhulst called this solution the "logistic curve ". The French term "lo- 
gistique " was used to signify the art of calculation. 

In a form of recursive equation it is named the logistic map [7] and was 
studied in [13], where Yorke has coined the term chaos. According to the au- 
thors, chaos is a result of presence of periodic orbits with all periods. Getting 
spectacular results the authors did not know about Sharkovsky's theorem |17j . 
Yakov Sinai [35] marks historical milestones of the chaos theory and anticipates 
its future: 

I believe that the future of the chaos theory will be connected with new phe- 
nomena in non-linear PDE and other infinite- dimensional dynamical systems, 
where we can encounter absolutely unexpected phenomena. 

This development existed long before the introduction of the term chaos. 
Vladimir Arnold [lj emphasizes the origins in the work of Henry Poincare. If 
we narrow the topic, then the orbits of a logistic map are given by H"6x(l — x). 
A chaotic time series is extremely sensitive to the initial conditions. Rounding 
off errors of computer operations evaluating initial conditions before a next 
iteration can become a natural source of instability in modeling such systems. 

Lorenz equations. Edward Lorenz [15] , using a modified model of Barry 
Saltzman |21j for simulation of atmospheric phenomenon, have discovered un- 
expected behavior of a solution of the system of three ordinary differential equa- 
tions referred today as the Lorenz equations 

^ = -oX + aY. i ^ = -XZ + rX-Y, ^=XY-bZ, 
dr dr dr 

where cr, r, and b are positive constants, and r is a dimensionless time. The solu- 
tion is a phase trajectory (X(t), Y(t), Z(t)), where X, Y, and Z define a phase 
space of a layer of fluid of uniform depth between two surfaces maintained at 
two different temperatures. Depending on the conditions and Rayleigh number 
the liquid remains steady or gets in motion, convection. X is proportional to 
the intensity of the convection, Y is proportional to the temperature difference 
between the ascending and descending currents, and Z is proportional to distor- 
tion of the vertical temperature profile from linearity. While the equations can 
be solved using different numerical methods, the author reviewed the double- 
approximation procedure described by Lorenz and found its iteral expression. 
Consider the initial value Pq of the vector variable P, and another variable p 



Po 



Then, a point of the approximated trajectory after the n iterations is equal to 
P n = W l po 0.5([ Yj+UUl °r -1 -z | \ I ) Ar+ ( y ))). 
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This is a direct way to a computer program. 

A financial time series. The most complete information about transactions 
comes from the markets as ticks, triplets of time, price, and number of contracts 
or shares named volume (tj,P,*,Vi). They are ordered by time. An ordinary 
daily session of a liquid nearby S&P500 E-mini futures contract traded elec- 
tronically on GLOBEX platform on the Chicago Mercantile Exchange brings 
several hundred thousands ticks. This associates with the name high-frequency 
trading [6\. These data is a result of cooperation of modern technology and 
human consciousness governed by partly unknown laws of nature. Substantial 
and frequent potential profits attract to trading. Real losses sober. Inability 
to avoid trading risk and the behavior of the market "attempting" to fool the 
majority of its participants creates a big assortment of approaches for "break- 
ing the market code" spanning from science to astrology. The markets appear 
random. A hope of traders is to find something more than a pure luck but a 
dependence between past and current behavior. Since all cannot win together, 
because the market does not create treasure but redistributes it giving away 
a portion to the industry supporting trading, non-zero sum game, the task of 
systematic winning by minor percent of participants becomes extremely tough. 

Probability theory, statistics, and theory of stochastic processes are tra- 
ditional ways to study the markets. Behavioral finance attempts to uncover 
psychological contribution into the process and find suitable quantitative meas- 
ures of it |30| . Forgetting about human being component of the activity can be 
a costly trading experiment. Albert Shiryaev calls the first order task study- 
ing statistics of waiting times, At, = i; — together with price increments, 
APi = Pi — Pi-i [24) p. 379 of Russian edition]. More and more empirical obser- 
vations confirm that Bachelier's or Samuelson's representations of price incre- 
ments or logarithms of price ratios, asset returns, by random Gaussian variables 
are only simplifications [20J. Robert Engle describes [5] that the Weibull [3"2"] 
distribution well approximates empirical waiting times. This as well as the gen- 
eralized Kumaraswamy distribution |12) are suitable for the waiting times com- 
ing from the Chicago Board of Trade [30] • Discrete probability distributions and 
laws such as multinomial, Hurwitz Zeta, and Zipf-Mandelbrot can approximate 
discrete price increments and their extreme values [20J. Prices and their incre- 
ments are essentially discrete and can be expressed as whole numbers of ticks. 
The author names the indecomposable further waiting times and price changes 
coming from a single trading session as a- and b-increments and the changes 
between the first, in a next session, and the last, in a current session, prices 
c-increments. These are quantitative measures of the a-, b-, and c-properties. 
Thus, the a-b-c-process evolves as following: a) the first property determines 
the moment of the next transaction in a trading session; b) the second property 
carries on the first one and determines the price fluctuation expressed by a whole 
number of ticks; and c) the third property is responsible for the discrete price 
change between the current last and next first prices. Selecting time Ut and price 
Up units for depicting an a-b-c-process on a chart, each tick becomes a point 
(tiUt,PiU p ). If this is viewed as a vector or a complex number with imaginary 
price component, then Azi = ((tt - U-i)ut,(pi — Pi-i)up) = {u t AU,u p Api). 
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Let Zo be the last tick of a previous session, AZ(AA, AC) be the random c- 
increment between a next first and a last previous ticks, and the lower case z 
be responsible for the intraday ticks. Then, using the iterals the a-b-c process 
is 

w J Zr H z ; AZ (z+Az), 

where j = 0, 1, 2, . . . enumerates sessions, Nj is a random number of intraday 
ticks, and Az is a random variable. To which degree everything is random, in- 
dependent, or varying in time is the subject of numerous modern investigations. 

A division of No on subsequences. A sieve. Below are the sequences 
of the nonnegative integers n, their positions p, and even e or odd o types t 

n 1 2 3 4 5 6 7 8 9 10 11 12 13 14 ... p ... 
p 1 2 3 4 5 6 7 8 9 10 11 12 13 14 ... p ... 
teoeoeoeoeo e o e o e ... e V o ... 

A sieve is applied to the sequence n. The sieve extracts numbers occupying 
even positions and puts them into the subsequence E, where their positions 
are renumbered in the natural order starting from zero. The extended types 
indicate even-even ee and odd-even oe even numbers located on even and odd 
positions within E. The E contains all and only even numbers including zero. 

n 2 4 6 8 10 12 14 16 ... 2p ... 

p 1 2 3 4 5 6 7 8 ... p ... 

t ee oe ee oe ee oe ee oe ee ... ee V oe ... 

The remaining numbers occupying odd positions in the original sequence n are 
put into the subsequence O, where they get new positions counted also from 
zero. The O contains all and only odd numbers. The odd numbers on even and 
odd positions within O get even-odd eo and odd-odd oo extended types. 

n 1 3 5 7 9 11 13 15 17 ... 2p+l ... 
p 1 2 3 4 5 6 7 8 ... p ... 

t eo oo eo oo eo oo eo oo eo ... eo V oo ... 

Being applied to E the sieve splits it into EE and OE. The new letter E or 
O, depending on the evenness of the extracted positions, is concatenated to 
a previous name on the left side. Similarly, from O the sieve produces EO 
and OO subsequences. The sieve is not applied to the subsequences getting 
the prefix EO in their names. These subsequences are referred to as the fixed 
subsequences. It is applied once to all other subsequences creating the following 
growing binary tree with stopping, final terminal, nodes corresponding to the 
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fixed subsequences. 



■■■■<- EE <- E <- N 

£0££ <- OEE EOE <r- OE EO <- O 

EOOEE <- OOEE EOOE <- OOE £00 <- OO 

EOOOEE <- OOOEE EOOOE <- OOOE EOOO <- OOO 



To save space the right child nodes are shown as down nodes. The three dots at 
the bottom indicate nodes, subsequences, to which the sieve is applied further. 
The tree dots on the top left indicate EEE to which the sieve is applied. 

After each application of the sieve to a parent's subsequence each member, 
number, is not missed and gets a new position in the left or right child's sub- 
sequences. Thus, after each finite number of iterative applications of the sieve 
every number is classified and occurs in one and only one terminal subsequence. 
Since 1) the sieve is applied to every E m E subsequence creating its right child 
node OE m E from all numbers on odd positions, 2) the position of any num- 
ber, except zero, alternates between even and odd after each new placement, 
and 3) the left nodes EOkO and EOE m E are created from each odd OkO and 
even OE m E subsequence, any number, except zero, after a finite number of 
iterative applications of the sieve becomes a member of a fixed subsequence. 
The numbers within a fixed subsequence are equivalent under the iterative sieve 
application. 

Weather a subsequence is fixed or not its name is unique. It also associates 
with a unique formula of infinite arithmetic progression linearly expressing any 
number within the sequence from its position counted from zero. Given a num- 
ber, a reverse linear function returns its position. The direct linear functions of 
p are in the nodes of the tree one-to-one mapping the tree of names. 

■■■•<— 4p <— 2p <— n = p 

\. \. \. 

16p + 4 <- 8p + 4 8p + 2 <- Ap+2 Ap + 1 <- 2p + 1 

\- -i- -i- 

32p+12 <- 16p+12 16p + 6 <- 8p + 6 8p + 3 <- 4p + 3 

-i- -i- -i- 

64p + 28 <- 32p+ 28 32p+14 <- 16p+14 16p+7 <- 8p+7 



Replacing the position p in the formula of a parent's node with the even 2p or 
odd 2p + 1 position produces the formula for the left or right child nodes. 

Geometrically, the points on the lines specified by the formulas associated 
with the fixed subsequences and computed for non-negative integer positions p 
do not coincide. This is true for terminal sequences on each iteration. 
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Iterals of linear functions computed for 2p and 2p + 1 initial values and dif- 
ferent number of iterations together with their nesting properties handle com- 
plexity. In order to get the iteral expression for a node of the tree the following 
common procedure is available. It is based on already described iteral properties 
and the fact that in a tree a unique path connects any two nodes. Our paths 
begin from No, which has the simplest linear function p. Making a move to a 
left node, apply V\\p to a parent's function. Making a move to a right node, 
apply H^p+i to a parent's function. At the start up, as we have seen in the 
section on iteral expression, the two identities are obtained 

E, 2p = Hl p p <- N -+ yil p+lP = 2p+l,0 

The iterals are nested on the path. As an example, consider the iterals nesting 
for the move N -> O -> OO -> EOO 

HljAl^JAl^p = Hj p Hj p+1 (2p + 1) = Mj p (4p + 3) = 8 P + 3. 

We shall be most interested in the fixed subsequences. They are subdivided on 
the subsequences EOkO containing only odd numbers 

Vll p V[ k 2p+1 (2p+l) = 2 k + 2 p + 2 k + 1 -l, fc = 0,l, ... , 
and subsequences EOiE m E 

nl (hu . . . K P+ ) i m;> = 2™+<+ 2 P + 2™+ i ( 2 ' - 1), 

m = 0,l,2, J = 1,2,3, p = 0,l,2, ... . 

containing only even numbers. Zero always remains in a subsequence of even 
numbers with the growing name EE^E, I = 0, 2°° +2 p,p = 0. By construction 
of the tree, we proved that each number becomes a member of one and only one 
of these subsequences. Thus, these iterals define the one-to-one correspondence 
between positive odd integers a and pairs of nonnegative integers (k,p) and 
between positive even integers b and triplets of integers (m > 0,1 > l,p > 0). 
In other words, V odd a e N 3 unique k,p e N solving the equation 

a = 2 fe+2 p + 2 fc+1 -l. 

Similarly, V even b g N 3 unique m,p g N A I g N solving the equation 

b = 2 m+l+2 p + 2 m+1 (2 l -1). 

While 6 = 0, 1 = 0, p = resolves the last equation for any m, we keep zero in 
EEooE. Given a positive number n the following algorithm consisting of three 
if-branches computes its k and p or m, I and p. 

If n = 0, then 1 = 0, p = 0, and m is arbitrary, set m = oo. 
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If n is odd, then subtract one and divide the result by two. If the result is 
odd, then repeat subtraction of one and division by two until an even number or 
zero is obtained. This is the stopping number. Count the number of repetitions 
on the way down to the stopping number: k = counter k — 1. Divide the stopping 
number by two. The result is p. Example, n = 39 — > 19 — > 9 — > 4. We stop at an 
even number, or zero, after making three subtractions of one followed by division 
by two. Then, k = 3 - 1 = 2. The p = stoppmg 2 "" m6er = § = 2. Verifying: 
2 k+2 p + 2 fe+i _ l = 2 2+2 2 + 2 2+i _ i = 39. Example: n = 7^3^1^0; k = 

3-1 = 2; p=§=0. Verifying: 2 fe+2 p + 2 k+1 - 1 = 2 2+2 + 2 2+1 - 1 = 7. The 
numbers 7 and 39 are equivalent. They belong to one fixed subsequence EOOO 
given by 16p + 7. The number 5 corresponds to (k = 0,p = 1) and belongs to 
EO given by Ap + 1 . The numbers 5 and 7 are not equivalent since they belong 
to two different fixed classes. Of course, both are members of O given by 2p+l. 
However, O is not a fixed subsequence and a subject for application of the sieve. 

If n is even, then it can be divided m+ 1 times by two giving the odd number 
2 l+1 p + (2 l — 1). The m = counter m — 1. The odd number is unit, indicating 
that given n is the power of two = 2 counteTm = 2 m+1 , only if p = and Z = 1. 
This means that all powers of two are extracted each by one and only one fixed 
subsequence of even numbers EOE m E and become their initial elements with 
positions p = 0. Thus, all powers of two are pair wise not equivalent. They are 
distributed between the fixed subsequences of even numbers as 

2 = 2*6 EOE : 8p + 2, 4 = 2 2 e EOEE : 16p + 4, ... , 

2 m+i e EO E m E : 2 m+3 p + 2 m+1 , .... 

For n = 2 m+1 the unique triplet is (m = counter - m — 1,1 = l,p = 0). When 
the odd number is not the unit, then we treat it as in the second if-branch 
above for the odd numbers and compute the counter k- The I = counter k 
and the p = ^ PP ing number Remember that k = counte r k - 1. Example, 
n = 28^14^7; m = 2-1 = 1; 7^3^1^0; Z = 3;p=§ = 0. Verifying: 
2 m+i+2 p+2 m+i(2i-\) = 2 1+3+2 p+2 1+1 (2 3 -l) = 64p+28 = 28. 28 is the initial 
member of the even fixed subsequence EOOOEE. The latter may not contain 
a power of two. Example, n = 38 -> 19; m = 1 - 1 = 0; 19 -> 9 -> 4; / = 2;p = 
| = 2.Verifying: 2 m +'+ 2 p+2 m+1 (2 i -l) = 2°+ 2 + 2 2 ?+2 0+1 (2 2 -l) = 16p+6 = 38. 
Thus, 28 € EOOE. The latter may not have a power of two. Its initial member 
is 6. 28 and 38 are not equivalent. 

The unique mapping <-» (m = 0, 1 = 0,p = oo), odd o (fc,p), and even o 
(m, Z,p) and subdivision on equivalence fixed subsequences allows reordering N 

< EOiE m E -< < EOE m E ~< 

< EOiEE -< < EOEE -< 

< EOiE -< < EOE -< -< EO -< < EO k O 

or with iterals 
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■■■<K P (hU) *v<---< K p K p+ ^ 2 p < 

-< Hj p (2p + 1) ■< ■ ■ ■ ■< Y[\ p V\ k 2p+l (2p+l)<--- 

All formulas for the iterals above are already derived. 

The formulas for odd a, and even b can be combined into one formula valid 
for any n?N, where m = — 1, 0, 1, . . . , I = k + 1 = 1, 2, . . . , and p = 0, 1, . . . 

n = 2 m+l+2 p + 2' n+1 (2 l - 1), n o (m,l,p), m G ZAm > G N,p G N , 
n = 2 m+i( 2 ;+i ?3 + 2 ; _ ^ 

All natural odd numbers are produced for m = —1, all natural even numbers 
are generated with m G No, and all powers of two, a subset of even numbers, 
are obtained for m € No A( = 1 Ap = 0. The formula does not miss any natural 
number and assigns a unique triplet (m, l,p) for the latter. We also can get zero 
for I = A p = and arbitrary m. Uniqueness does not hold for zero. 

Waclaw Sierpinski calls Mersenne numbers the numbers represented as 2 l — 

I, 1 = 1, 2, ... [26"1 p. 76-79]. The subset of prime numbers, Mersenne primes, are 
named after the Minimite friar, Farther, Marin Mersenne |2] p. 334, 385], [El p. 

II, 22]. Computer specialists working with the standard binary representations 
of integers well know that 2 l - 1 = 1 x 2'" 1 + ... 1 x 2 2 + 1 x 2 1 + 1 x 2° = 
S^=o _1 ^ ■ For instance, decimal 7 is equal to binary 111 = 1 x 2 2 + l x 2 1 + 1 x 2°. 
We see that each fixed arithmetic progression for odd numbers 

2 fc+i_ 1 _ 2' +1 p+2'-l = 2 l+1 p+Mi starts from a Mersenne number Mi and has 
the progression difference equal to the power of two obtained after incrementing 
by one the position, starting from one, of a number within the sequence of 
Mersenne numbers. The formulas for the fixed arithmetic progressions of even 
numbers are obtained from the odd counterparts after multiplying by a power 
of two. 

All spectrum of interesting mathematical questions related to the infinite 
subsequences can be applied to the equivalence arithmetic progressions con- 
sidered in this section. For some of them, such as EO : 4p+l, the representation 
of the prime numbers of the sequence by the sum of two squares have become 
classics of the number theory |26l p. 51-53], [SJ p. 20]. How the prime and/or 
Mersenne numbers are distributed between the fixed subsequences? It is worth 
to quote Sierpenski's words (author's translation to English) [26l p. 80]: 

Question, whether a given infinite sequence, defined even in a simple way, 
contains infinitely many prime numbers, in general, is quite difficult. 

The task to show that iterals, iteral notation, and the sieve are natural means 
to derive the equivalence arithmetic progression formulas is completed. 

The subdivision of No is interesting for studying the division by two applied 
to the even numbers in the even fixed subsequences on the left side of zero in 
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conjunction with the operation 3a: 2 f 1 applied to the odd numbers in the fixed 
subsequences in the right side of zero. This is done in the next sections. 

The operation 3a: 2 l ~ 1 applied to EO^O. There is unsolved convergence 
problem named after Lothar Collatz Q3j. If the unit is given, then stop. If an 
even number greater than one is given, then divide it repeatedly by two until the 
unit or a greater odd number is obtained. If it is one, then stop. If the result 
is an odd number greater than one or an odd number is given initially, then 
apply to it repeatedly 3a ^ 1 " 1 . If an even number is obtained at any step, then 
switch to division by two. The unproved proposition is that for any positive 
integer after the finite number of steps the result will be the unit. For instance, 
3->5->8->4->2->l. 

Within the scope of this article it is natural to denote the result of 3x ^ 1 
applied k = 0, 1, . . . times by the iteral 3a: 2 1 ' 1 . In order to get integer, the 
operation must be applied to an odd number n — 2p + 1 . Under this condition, 
one iteration is equivalent to adding to an odd number a position of a next odd 
number within O. Indeed, 

Hi 3x + I 6v + 4 

The result is an odd number, if p is odd. Otherwise, the result is an even integer. 
We just proved that in terms of our subsequences, the operation converts an 
even-odd number to an even number: EO — > E. Getting an even number from 
an odd one is important for the entire proposition because this is the only way 
eventually to face with an even number, which is a power of two leading to final 
one. Unfortunately, the same operation converts an odd-odd number to an odd 
number. Let us prove that being applied k + 1 times to an odd number from the 
fixed subsequence of the order k EOkO an even number is obtained. Indeed, 
for the operation applied j times to such a number we get 

W , , 3x + l T/p 3x + l 

= 3 J '2 fc+2 - i p + 3 i 2 fc+1 -^ - 1. 

As long as j < k the result remains an odd number. When j achieves k + 1 
the result is even and equal to 3 k+1 2p + 3 k+1 - 1 = 3 fe+1 (2p + 1) - 1. Thus, 
any odd number n O {k,p) after k + 1 iterations applying the operation 3x + 1 is 
converted to an even number larger than n. It is larger because on each iteration 
it is added by increasing position of a next odd number. 

What are (k new ,p n ew) after application of 3:r 2 ^ 1 to an odd n «-» {k,p)l 



2 fe + 2 p+2 fe+l_l 



^±l = 3 1 2 fc + 1 p + 3 1 2 fc -l = 



= (2 + l)2 fe+1 p + (2 + l)2 fe - 1 = 2 k+2 p + 2 k+1 p + 2 k+1 + 2 k - 1. 

Let us apply the algorithm for odd numbers to find k new and p new . The right 
most —1 is not a problem because we subtract one and then divide by two until 
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even or zero stopping number. This will keep —1 at the end of expression intact. 
Thus, after counter = k steps the stopping number is even 2 2 p + 2p + 2 + 1 — 1 = 
2 2 p + 2p + 2. Then, k new = k — 1 and p new = 3p+l. This operation reduces the 
order of the odd fixed subsequence EOkO —> EOk-iO and increases position 
within the new subsequence: (k,p) — > (k — l,3p+ 1). After k iterations the 
converted number is still odd and belongs to EO : Ap + 1 with the order 0. It 
is one step from becoming even. On jth application of the operation on the 
way to smaller orders Vj < k the odd number is given by the pair (H^(fe — 1), 

H£(3p + 1)) = (k — j, Vp + YllZo 3*) = (fc - j, &P + ^r)- The sum for j = 
is set to zero as it is described in the historical section for J^. Here are the 
starting and ending (ki ast , piast) nodes on the iteration way, when the numbers 
are still odd 

(fc.p) G EO k O : 2 k+2 p + 2 k+1 - 1 -> (0, 3 k p + ^-y-^) €EO:4p+l. 

Single application of the operation strips off one O from the name of an odd fixed 
subsequence. The (fc + l)th application returns an even number 3 fe+1 (2p+l)-l. 
It can be obtained as we did it above or from 4p + 1 replacing p by pi as t and 
applying the operation the last (k + l)th time, when it still can be applied. 

The operation | applied to EOiE m E. The even numbers (m,l,p) from 
these subsequences are given by 2 m+l+2 p + 2 m+1 (2 z — 1) and can be divided 
(m + 1) times by two 

FT +1 x -2'+ 1 « + i 

The result is an odd number (k,p) from a subsequence EOkO. Indeed, for 
I = k + lwe get 2 k+2 p + 2 k+1 — 1. Being applied (m+1) times to a number from 
an even fixed subsequence, division by two strips off all suffixes E: EOiE m E — > 
i?0; = EOi-iO = EO^O. The original even number and the obtained odd 
number have the same numerical position each within its own subsequence. The 
intrinsic order of an odd number "hidden" in the structure of an even number 
does not change during this transformation. With respect to division by two 
I = k + 1 and p are invariant. 

Using the current layout of the tree diagram, where the right child nodes 
are shown as the bottom nodes, division by two starting from an even fixed 
subsequence is a horizontal move to the right, sequentially visiting nodes of 
other fixed even subsequences with incrementally reduced number of E in the 
suffix of the name. This continues until the node of an odd fixed subsequence 
with the same number of Os is visited. Then control is given to 3x + 1 , which in k 
vertical moves visits EO. The latter is a subject of the last (fc + l)th application 
giving an even 3 k+1 2p + 3 k+1 - 1. Which even subsequence EOi ni!m E mni!m E : 
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2 m ncw +i ncll ,+2 pnew + 2 m„ e .+i( 2 l» e » _ i) does it belong to? 

EO 

i/l, ^r 1 

EOtE-iE = EO lnew E mnew E ' t fc = i - 1, 

EO k OE m E -> £O fc O 

m + 1, divisions by 2 

We have proved earlier that 2 m+1 G EOE m E : 2 m+3 p + 2 m+1 is the initial 
p = member of the subsequence. In order the main conjecture would be true, 
eventually, a move from EO should visit one of these subsequences and exactly 
its initial member: l new — lAp new = 0. "Eventually" emphasizes that this does 
not have to happen on the first loop. Then, both endless 1) cycling between 
EO and non-initial member of EOE m E and 2) general loops must be excluded. 

The fate of even number 3 fe+1 2p + 3 k+1 — 1. The even if-branch of 
the algorithm determining (m,l,p) begins from division by two until an odd 
number is obtained. If the odd is one, then the source is the power of two, 
otherwise, m = counter m — 1 and the odd if-branch of the algorithm takes 
control. The k and p are arbitrary and independent because any number can be 
given at the beginning. Without loosing generality, there are four possibilities: 
1) fc = A p = 0, then the even number is two and we do not have to be here 
because the original given number is 1 £ EO : 4p + 1, indeed, 2 = I ; 2) 
k = OAp = odd p 2 q , oddp > 1, q > 0, then the number is 6p + 2 = 3odd p 2 q+1 +2 
and Vq > m new = 0; 3) k = odd^2 r hp = 0, odd^ > 1, r > 0, then the number 
is 3 fc+1 — 1; and 4) k = oddk2 r A p = odd p 2 q producing the most complex 
expression, where we apply the Newton's binom |10) for the second occurrence 
of 3 fe+1 = (2 + l) fc+1 . Euler's an d Kramp's ! are in use 

3 k+1 2p + (2 + l) fc+1 - 1 = 3 k+1 odd p 2 q+1 + Y ( fc + 1 ) ! 2 fc+1 ^ - 1 = 

3 odd v 2 +^ jl{k + 1 _ jy 2 
1 k+1 ndd 2 q+1 4- °V {odd k 2 r + l)l oddk2 r +1 ^j 

p + jKodd^ + i-jy 2 

The (k + l)th term 1 is taken out of the and the last value of the index j 
is reduced to k. The (fc)th term of the sum is equal to (k + 1)2 = (oddk2 r + 
1)2 1 . Comparing with other terms, for r > it contains the smallest factor 
represented by a positive power of two. For r > the sum is divisible by 
two only once. By the same reasoning in the third case for r > the number 
3 +1 — 1 is divisible by 2 only one time. Thus, \/odd p , oddk, {,r£N m new — 
and3 fe+1 2p + 3 fc+1 -l G EOi new E = EO hnew+x E : 2 l —- +2 p new + 2{2 l ^~l) = 
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2 kncw+3 p new + 2(2 kncw+1 — 1). After single division by two we get for the second, 
third, and fourth cases 

3odd p 2 q + 1, odd for q > and even for q = 0; 

oddu2' 



y {odd k 2 r + iy. ddkr _ odd and eyen for = 

^ j!(odd fe 2- + l-j)! 



^ 3 } -(odd k 2 r + 1 - j)! 

In the fourth case the number is odd for (q > A r > 0) V (q = A r = 0) 
and even for (q > A r = 0) V (q = A r > 0). Determination of m new is 
not over for even result. While for any given natural number n the algorithm 
returns (m > 0, 1 > \,p > 0) in a straight forward computation, without an 
explicit expression for EOi new E mnew E : 2 m ^ +l ^ +2 p new + 2 m — + 1 (2 / "=™ - 1) 
it is difficult to move forward. The conjecture remains unproved. It is proved 
that any natural number, which is not a power of two, will be converted to an 
even number using the finite number of divisions by two and 3a: 2 t " 1 operations, 
which can be determined in advance by setting the correspondence between the 
given number and pair (k,p) or triplet {m,l,p). Application of iterals to these 
operations looks natural. 

The author is grateful to Florian Cajori, who many years ago has written 
the words suitable for ending this article |3, p. 77] 

We know that intellectual food is sometimes more easily digested, if not taken 
in the most condensed form. It will be asked, To what extent can specialized 
notations be adopted with profit? To this question we reply, only experience 
can tell. It is one of the functions of the history of mathematics to record such 
experiences. 



Acknowledgments. I would like to thank Alexander Tumanov and a referee 
for reviewing the article and useful advices and criticism. Due to them I have 
added the second version of the iteral notation, where in ambiguous cases a 
variable of a function is specified explicitly. 



6 Appendix 

The author has written a program oneness. The CH — h code can be requested by 
email. Given a positive integer and radix of the system (decimal, octal, binary, 
hexadecimal, and, in general, from two to 36) it returns intermediate numbers 
for a Collatz's task. The program works with the numbers not exceeding in 
computations the Mersenne number 2 31 — 1 = 2147483647. Here is the output 
for the run: oneness 9 3 
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4 : (m=l, 1=1, p=l) 




20 


8 


10 


D2 


101 


E0E : 8p + 2 


: (m=0, 1=1, p=l) : 


10 




9 


5 


D2 


12 


E0 : 4p + 1 : 


(k=0, p=l) : 5 






10 


8 


3X 


22 


E0E2E : 32p + 


8 : (m=2, 1=1, p=0) 




8 


11 


4 


D2 


11 


E0E1E : 16p + 


4 : (m=l, 1=1, p=0) 




4 


12 


2 


D2 


2 


E0E : 8p + 2 


: (m=0, 1=1, p=0) : 


2 




13 


1 


D2 


1 


E0 : 4p + 1 : 


(k=0, p=0) : 1 







The first column is the step counted from zero. The second column is an inter- 
mediate number, where nine on the line zero is the input. The third column is 
the name of the operation, where IN is input, and 3X or D2 is £ or division 
by two applied to a previous number. The fourth column is the number repres- 
ented with the radix three requested in the command line. The last descriptor 
with the colon delimiters contains the name of the subsequence, the formula of 
the arithmetic progression, the pair (k,p) or triplet (m,l,p), and repetition of 
the same number computed after reusing the pair or triplet to make sure that 
the program works correctly. 
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